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The recently developed formalism of Markovian master equations for quantum open systems 
with external periodic driving is applied to the theory of dynamical decoupling by periodic control. 
This new approach provides a more detailed quantitative picture of this phenomenon applicable 
both to traditional NMR spectroscopy and to recent attempts of decoherence reduction in quantum 
devices. For popular choices of bath spectral densities and a two-level system, exact formulas for 
the relaxation rates are available. In the context of NMR imaging, the obtained formulas allow the 
extraction of a new parameter describing bath relaxation time from the experimental data. 
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I. INTRODUCTION 


The historical origin of dynamical decoupling dates back to the beginning of 1950-ties when the “spin-echo” or 
“refocusing” phenomena were observed for nuclear magnetic resonance (NMR) [1]. Namely, a properly designed 
sequence of pulses, which mimicks the time-reversal operation, reduced the decay of coherence for the NMR signal 
emitted by a large sample of nuclear spins. The “static” source of decoherence due to the inhomogeneity of the 
external magnetic field dominated over the “dynamical” decoherence caused by the interaction of a single spin with 
environmental noise. Therefore, the effective elimination of the former allowed the high-precision measurement of 
decoherence time T 2 due to the latter. Further progress in experimental techniques allowed the shortening of the 
repetition time of the pulses T below the values of T 2 , leading to the reduction of dynamical decoherence as well [2]. 
In Ref. [2] the elementary theory of this effect, called coherent averaging , has been developed. 

In the last two decades new advances in the fields of quantum information, quantum computing or, more generally, 
quantum control have led to the revival of this topic both on the theoretical and experimental sides (see Ref. [3] 
and the book [4]). Recent theoretical contributions concentrated mostly on the development of numerous, often quite 
sophisticated, decoupling schemes such as group-based designs, Eulerian design, or concatenated schemes which went 
far beyond simple ’’‘bang-bang’” periodic decoupling. However, an approach based on coherent averaging applied to 
a generic bath could not give detailed quantitative predictions concerning the decay rates except for the universal 
condition for decoherence suppression by periodic control, 
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where cj cut is supposed to be the cut-off frequency for the spectral density of the bath. Much more detailed information 
has been obtained for the exactly solvable case of Gaussian classical noise combined with ir- rotations and for various 
decoupling schemes; see, e.g., Refs. [5-7]. 

The inequality (1) was often treated as a sign of the non-Markovian character of the setting. Another picture 
invoked in the context of dynamical decoupling was the analogy to the quantum Zeno effect- the freezing of dynamics 
due to frequent measurements. Perhaps the most insightful, from the physical point of view, is an explanation based 
on the Fermi golden rule. Here, a frequent periodic perturbation shifts the Bohr frequencies of the system into a 
frequency domain higher than that corresponding to relevant excitations in the bath [8]. As a consequence, energy 
conservation suppresses efficient coupling to the bath. 

The aim of this paper is to apply the recently developed theory of quantum Markovian master equations for 
periodically driven systems to dynamical decoupling by periodic control operations. This framework introduced in 
Ref. [9] and studied in Ref. [10] was applied to a number of models of quantum engines and cooling schemes [11-15]. 
It is based on the combination of Davies’ weak-coupling-limit technique [16] with Floquet’s theory of periodic-in-time 
quantum Hamiltonians [17]. The main result of the formalism is the master equation written in the following form: 
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where pi n t(t) is the time-dependent reduced density matrix of the system in the interaction picture and C generates 
a completely positive dynamical semigroup, i.e., possesses the canonical Lindblad-Gorini-Kossakowski-Sudarshan 
structure [18, 19]. The latter property means that the dynamics is Markovian, which seems to contradict the popular 
opinion about the non-Markovian nature of the dynamical decoupling setting. 

The origin of this controversy is again a very popular picture of the Markovian approximation invoking the short 
correlation time of the bath as a necessary condition. Actually, this correlation time can even be non-well-defined, 
as in the case of spontaneous emission where the electromagnetic bath correlations decay like t~ 4 but, nevertheless, 
spontaneous emission is a perfect example of Markovian (exponential) decay. A more rigorous analysis shows that 
Markovian behavior is a consequence of the following conditions: (i) weak coupling to a large stationary bath with a 
quasicontinuous spectrum and approximately Gaussian correlations, (ii) a mild assumption like integrability of two- 
point correlation functions of the bath implying a continuous spectral density function, slowly varying on the scale 
of typical relaxation rate, and (iii) well-separated relevant Bohr frequencies of the system leading to an efficient time 
averaging of the Hamiltonian-induced oscillations during the typical relaxation time. 

As a consequence of the above conditions one should treat the equation (2) as describing the coarse-grained-in-time 
dynamics valid for times longer than the timescale determined by the Hamiltonian dynamics. In the case of fast 
periodic control this timescale is determined by the pulse repetition time. A rule of thumb for the applicability of the 
Markovian approximation is an essentially continuous energy spectrum of the bath and relaxation rates essentially 
lower than the relevant Bohr frequencies of the system, including pulse frequencies. 

The main advantage of the proposed formalism is the existence of explicit formulas for the effective decoherence 
and relaxation rates for periodically controlled systems depending on the details of the free-system Hamiltonian, the 
form of perturbation, the system-bath interaction, and the parameters of the bath. In simple but important cases 
even closed exact solutions can be obtained that allow a rigorous analysis of dynamical decoupling phenomena. 


II. OPEN SYSTEM WITH PERIODIC DRIVING 


In this section we briefly present a theoretical framework of periodically driven open quantum systems, first in 
a general case (in Sect. II A), then in a case of periodic kicking (in Sect. II C). We skip most of mathematical 
technicalities; the reader can find more details in Refs. [9, 10, 15]. We also use the same units for energy and angular 
frequency, i.e., h= 1. 

We consider an open quantum-mechanical system S described, for simplicity, by a finite-dimensional Hilbert space 
and by a time-periodic Hamiltonian H(t) = H(t + T). At this point we note that H{t) must be understood as a 
properly renormalized, physical Hamiltonian, already containing all Lamb-like shifts induced by interaction with the 
environment [10, 20, 21]. The system-bath interaction is parametrized as 

Hint = ® Rg, ( 3 ) 

a 


where S a and R a are Hermitian operators of the system and bath, respectively. The initial state of the bath (with 
the average denoted by (•••)#) is invariant with respect to the free dynamics of the bath and satisfies (. R a )s = 0. 

Applying the standard derivation of the master equation in the interaction picture based on the weak-coupling-limit 
technique, one obtains equation (2) with the generator 

Cp = EE 7 a/j(w') (s a [w')pSp{^ - (</),/>}) • ( 4 ) 

fw'} ' ' 


Here, operators S a (u f ) are defined by the Fourier decomposition 


where 


s a (t) = u(t)\t)s a u(t) = 


jU{t) = —iH(t)U(t), U(0) = I, 


( 5 ) 

(6) 


and a set of frequencies {uj'} will be determined later using Floquet theory. The relaxation and decoherence rates are 
determined by the Fourier transforms of bath correlations which, for any c j, form a positively-defined matrix. 


oo 

7a/s(y = J e.~ lUt {Ra(t)Rp). B dt. 

— oo 
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Notice, that in both expressions (5), (7) and in the similar ones later on, the time dependence of observables is always 
given in the Heisenberg picture, governed by internal Hamiltonian of the system or bath, respectively. The formula 
(7) shows also why a short correlation time for the bath is not a necessary condition for the validity of the Markovian 
approximation. What is important is that the integral restricted to the interval [—r, r] converges to the definite value, 
describing the relaxation matrix element. This convergence can be fast even for slowly decaying correlations if only 
the relevant Bohr frequencies {cj} are high enough. 


A. Application of Floquet theory 

Applying Floquet theory to the unitary propagator [/(£), one obtains the following decomposition [10]: 

u(t) = P(t)e~ iSt , (8) 

with P(t) being periodic and unitary; H is a Hermitian averaged Hamiltonian , satisfying the relation 

U(T) = e~ inT . (9) 

The Floquet operator U (T) is a unitary operator, which propagates a quantum-mechanical state 2 p(t) over the full cycle 
[0, T\. Both U(T) and the averaged Hamiltonian H posses a common eigenbasis {</>&}; therefore, U(T)(j>k = e~^ ekT (j)k, 
where } are called quasi-energies of a system with a periodic Hamiltonian. This structure implies a particular form 
of Fourier decomposition (5) given by 


S a (t) = E E q)e i{u+qQ)t (10) 

q£Z {a;} 

where operators S a (uj,q) are subject to the relations 

5'a(w,?) t = S a (-LJ t ~q), 

[H, S a (u>, q)] = u>S a (u>, q). (11) 

The set of frequencies {c</} in Eq. (5) is labeled by two elements {c</} = {uj + qd; q — 0, ±1, ±2, ... } where 
{cj} = {e/e — e/} are Bohr-Floquet quasifrequencies, and Q = 2tt/T . The harmonics qCl correspond to the energy 
quanta qQ which are exchanged with the external source of periodic driving. 


B. Properties of Master equations for periodic driving 

By using the decomposition (5) one can rewrite the interaction-picture generator (4) in a the form 

£p = EEEo'«/3( w + 9^) { S a (uj,q)pSp(uj,q ) f - ^{Sp(u>, q) f S a (u, q), p} j . (12) 

a,(3 qEZ {u;} ^ ' 

From the properties (11) it follows that the generator (4) independently transforms diagonal and off-diagonal elements 
of p, computed in the eigenbasis of the averaged Hamiltonian H. 

It is often useful to employ the Markovian master equation in the Schrodinger picture which possesses the following 
structure: 


^ = -iim, m + (w(t)£ w(t)t) P (t) (is) 

with the unitary superoperator U(t)p = U(t)pU(t)^. The properties (11) imply also a very useful factorization 
property of the dynamics governed by (13) 


p(t) =U(t)e tC p( 0), 


(14) 


which allows to discuss separately the decoherence or dissipation effects described by C and the unitary evolution 
U{t). 
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C. Periodic kicking 

Now we focus on a very special case of periodic Hamiltonian; namely, the so-called kicked Hamiltonian , which can 
be formally given as 


H{t) = H 0 + \W 5 ( t ~ fcT )’ (15) 

k= — oo 

where A is a cumulated kicking strength and Hq and W are Hermitian operators. As was already suggested, the key 
point for open-system analysis lays in the Floquet operator U(T) and its spectrum, which allows one to write the 
appropriate formula for the interaction-picture generator C (4). A formal expression for the Floquet operator in the 
case of periodically kicked systems has been known for quite a long time already (see, e.g., Refs. [22-25] and references 
therein). Because of the singular nature of H(t), heuristically speaking, one is interested in a unitary evolution taking 
the state 'i/jft) not from time 0 to T, but rather from a time 0 + just after the time of the first kick to a time T + 
just after the time of the second kick. In the infinitesimally short interval [T _ ,T + ], the system experiences a kick 
of infinite strength by the operator AIF, so, effectively, the influence from Hq can be omitted. As a consequence the 
Floquet operator is given by 

U(T) = e ~ iXW e ~ iH o T = e ~ iRT . (16) 

The time-dependent unitary U ( t ) can be expressed in terms of two functions {x} and [x\ , called the fractional part 
and the floor function, respectively (see the Appendix), i.e., 

[x\ = max {k E Z : k ^ x }, (17a) 

{x} = x — [x\ . (17b) 

The {x} function is also known as the sawtooth function and is periodic with period 1, which directly implies that 
function {|r} will be of period T. But this allows one to write, 

U (£) = e ~ lHQT { ? } e^ T { ^ } e~ l ^ t (18) 

where H is given by (16). This is naturally consistent with the representation (8) implied by the Floquet theorem, such 

that P(t) = e~ lH ° T { ^ } e ll * T { t } is a periodic time-dependent operator; because the sawtooth function is periodic, 
this guarantees that results obtained in Ref. [10] apply and the whole construction presented in Sec. (II A) can be 
employed. It is important that, from the very construction, function U(t) given as Eq. (18) is discontinues at times 
t = nT, which is a direct consequence of the singular form of the Hamiltonian. Therefore, the resulting completely 
positive dynamics Aft) = U(t)e tC will be discontinuous as well. 


III. EXAMPLES OF ’’BANG-BANG” CONTROL 

In order to illustrate the power of the above formalism we discuss the simplest, but practically very useful, model of 
a two-level system (TLS) which can be physically realized as, e.g., spin 1/2, superconducting qubit, quantum dot or 
fluorescence center in solid. We restrict ourselves to the most popular case of periodic control by short pulses which 
can be treated within the formalism of kicked dynamics. The chosen details of the Hamiltonian and the interpretation 
of the model parameters are motivated by the NMR applications but obviously can be extended to other examples. 

We consider a two-level system described by the standard Pauli matrices. The system unperturbed Hamiltonian is 
given by and the control is executed by the external field oscillating with the angular frequency cj ex t and with 

the envelope /(£). The total TLS Hamiltonian in the rotating wave approximation reads 

H'ft) = ^cjocr 3 + fft) (e~ iuJextt (j + + e iuJextt cr~) . (19) 

The generic interaction with the environment can be written as 

Hint =v-R (20) 

where R = R 2 , R 3 ) represents three observables of the environment which in the NMR context, are usually chosen 

to be statistically independent. Then each component of the interaction can be treated independently and the cross 
correlations between different Rj can be neglected. 
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A. Example 1: Two-level system with “longitudinal” coupling 

First, we consider the system-bath interaction Hamiltonian (3) of the form 

H int = a 3 ® R 3 . (21) 

This type of coupling, called in the NMR context the spin-spin interaction , often dominates and hence, is interesting 
to be studied separately from the others. 

In the absence of external control, i.e., for f(t) =0, this model describes the pure decoherence of the TLS with the 
decoherence time traditionally denoted by X \ and given in the Markovian approximation by 

oo 

= 7||M= I e- iut (R 3 (t)R 3 ) B dt. (22) 


In order to simplify the analysis we describe the system in the “rotating frame,” i.e., we perform a unitary trans¬ 
formation Uq (t) defined as 


Ul(t)A = e hit^<r 3 Ae -±itu, ext * 3 ( 23 ) 

for an operator A, which will effectively switch us to the Heisenberg picture with respect to ^co ex t (7 3 . In this rotating 
frame the dynamics will be described by transformed Hamiltonian 

H"{t) = ^cr 3 + f{t)a\ (24) 

where 

A LOq ^ext (25) 


stands for the detuning parameter. 

Now we assume that the envelope f{t) is an infinite sequence of short pulses of duration much shorter than their 
separation time T which is, on the other hand, much shorter than the decoherence time XE Therefore, we can replace 
pulses by Dirac deltas to obtain the final effective TLS Hamiltonian in the rotating frame 


H(t) = ^a 3 + A 


E S(t-kT) 


,k= — 00 


a . 


(26) 


We choose A = tt/2 (“magic angle,” 7r-rotation) what drastically simplifies the description, but the exact solutions, 
albeit very involved, are also available for any A. The Floquet operator U(T) is, according to Eq. (16), given as 


and the Floquet basis is 


U(T) = 




iT A , 
2 C 


0 

0 — hi AT 


e2 


iAT 


—iTH 




1 

—es 


iAT 


1 


(27) 

(28) 


(29) 


Quasienergies, i.e., eigenvalues of iL, are found to be ei ? 2 = ±7r/(2X) and quasifrequencies are simply uj± = ± 7 r/T. 
The unitary propagator U(t) can be easily found by direct application of Eq. (A6); in the Floquet basis {^ 1 , 2 } it is 
given by 


/ e M#Jcos(^{£}) ie-MfUsin^q,}) \ 

( ie? i7T \-T\ sin (4p {£}) cos (4p {£}) ) 


U(t) = 


(30) 
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This allows to obtain the Heisenberg picture evolution of S' = cr 3 , 

oo 

<r 3 (f) = J2 ( y S(+;q)e i (T+^)t + s ^_.q^ e -it(f+ q n)t^ 
q= — oo 

with matrices 

* 5 (+; q) = 

S(-,q) 


2 i 


0 0 


7r(2 q + 1) \ 1 0 

2 i /0 1 


7r(2 q — 1) \ 0 0 

We can cast Markovian master equation (2) into an infinite sum 

oo 

£p = ^ ' £-qP 


(31) 

(32a) 

(32b) 

(33) 


q =—oo 


where C q is defined as 


( 2 7((j+g)^) . 2 7((~g+g)^) 


£qP — o 


(l+2q) 2 


Pll + 


(l-2q) 2 




_ ( 7((~i+g) Q ) , 7((i+g)n)\ 

l ( l - 2 ?) 2 + ( 1+ 2 ?) 2 J 


7(-(J+g) fi ) , 7((J+«) fi )W \ 

" p 22 l (l-2q) 2 + (l+2q) 2 J P 12 

2 7((l+g) fi ) 2 'i((~h+ < i) n ) 

P 21 (l +2g ) 2 PH - (l-2 g ) 2 W J 


(34) 


The form (33) raises a natural question regarding summability of the infinite series; however it is easy to show that 
since y(w) ^ 0 the series indeed converges whenever j(oj) is bounded. 

We consider now the Lorentzian spectral density corresponding to the exponential decay of bath correlations with 
the decay time r c 


7||M = 


1 


T 2 1 + T 2 UJ 2 ’ 


(35) 


where T 2 is the standard pure decoherence or dephasing time measured in the absence of modulation. The formula 
(35) is a reasonable high-temperature approximation for the “spin-spin” interaction with the bath, often used for 
NMR. 

With the choice (35) the summation in the formula (33) can be executed and produces the interaction picture 
generator 


£p= ~V\\( T ) 


pll — p 22 P12 

P 21 P 22 — Pll 


where the relaxation rate p\\{T) is expressed by the formula 


^^(i-^tanhT). 


(36) 


(37) 


Figure 1 shows a plot of the relaxation rate r]\\ as a function of kicking frequency Q compared to the Lorentzian 
spectral density (35). 

Now one can solve the interaction picture Markovian master equation and then transform the solution p int (£) = 
e tC p{ 0) back to the Schrodinger picture in the laboratory frame to obtain the exact expression for system’s density 
operator in the original eigenbasis of a 3 operator, 


Pit) =U 0 (t) (U(t)e tc p( 0)) 

with explicit matrix elements pij(t ), i, j G {1,2} given by formulas (77 = 77||(T)) 

pu(t)= 1 - [l + (-l)We-^( Pll (0)-p22(0))] , 

P2i(t) = ^e- 2 ^e-i( AT L*J““oO (e- iAT p 12 (0) +p 21 (0)) 

- i (— 1 ) LtJ e _, ? t e _ *( AT (i+ LtJ )-‘ J o t ) (p 12 ( 0 ) - e iAT p2i(0)) 


(38) 

(39a) 

(39b) 


P22 (t) = 1 - pn (t ), P12 (t) = p2 1 (f) • 


(39c) 
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V 0 ) [s ' 1] 



FIG. 1. Comparison between the relaxation rate rj y(Q) (black) and spectral density 7 y(Q) (gray, dashed) for the longitudinal 
coupling. 


We rewrite the solution (39) in terms of Bloch ball parametrization for the TLS density matrix 


p=^[I + x-a\ (40) 

with components Xi{t) of the Bloch vector given explicitly as 

xi = pu + P 21 = 2Rep 2 i, (41a) 

X 2 = i(pi 2 ~ p 2 i) = 2Imp 2 i, (41b) 

X 3 = Pu - P 22 = 2pii - 1. (41c) 

By using the definition 

4>(t) = w ext f + AT |^| , (42) 

where {x} c = {x} — 1/2 can be called the centered sawtooth function , one obtains the solution of Eq. (39) in the form 

x\(t) = e~ 2r]t cos0(£) (xi(0) cos 0(0) + # 2 ( 0 ) sin 0(0)) (43a) 

+ ( — 1) LtJ e _??t sin (j)(t){x i(0) sin0(0) — £ 2 ( 0 ) cos 0(0)) 

X 2 (t) = e~ 2r]t sin0(t) (^ 2 ( 0 ) sin0(O) + £i(0) cos0(0)) (43b) 

+ ( — 1) LtJ cos cj)(t){x 2 ( 0 ) cos 0(0) — £i(0) sin 0(0)) 
x 3 (t) = (- l)7J e -^ 3 (0). (43c) 


One can decompose the evolution of Bloch vector into three types of motion: (a) rotation with the angular frequency 
c^ext perturbed by the phase modulation of a sawtooth shape, (b) exponential decay with two decay rates, and 

277|| (T), an d (c) inversion of the slowly decaying component performed at the times nT. The exponential decay of all 
components of Bloch vector drives the TLS to the final maximally mixed state. 

1 . Spin-echo phenomenon 


In real experiments one measures the NMR signal produced by the x\{t),X 2 (t) components averaged over a large 
sample of individual nuclear spins. Because the external magnetic field which defines the Larmor frequency ujq is 
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not perfectly homogeneous the formula for phase (42) contains a deterministic component cj ext £ and the essentially 
random detuning A = c^o — cj ex t- Averaging the terms, which describe rotation in the X 1 X 2 - plane, with respect to 
fluctuations of A one obtains 


(cos </>(t)) = (cos (at jtj ) ) cosoJexti - ^sin (AT j t j ^sinwexti 


(sin <j>(t)) = (cos (AT 11) j ^ sinw ext £ 


+ (si "H4) 


COS CJgxt^* 


For times t = (n + 1/2)T the function {^} c is equal to zero, which implies that 


cos (j) | ( n + TJ \ = cos fcj ext (n+^)T 


and similarly 


(44) 


(45) 


(46) 


sin ( 


U ^i) T ) ) = Sil1 \ GXt \ ^ ^ T ^ 


(47) 


Hence, the effect of magnetic field inhomogeneity is undone in those particular moments of time. This is the spin-echo 
phenomenon and the formulas (42)-(45) allow us to compute the shape of the NMR signal for a given model of 
magnetic-field variations. 


2. Freezing of dynamics for fast pulses 

Due to the spin-echo phenomenon, in NMR experiments, for the pulse period T t c , the measured decoherence 
time is equal to T^(T) = rj\\ (T) _1 ~ T 2 , where X 2 is the ideal decoherence time for a single spin in the absence of pulses 
[see formula (37)]. Increasing the frequency of pulses to the values higher than r~ x (T r c ), one observes suppression 
of the relaxation rate r}\\(T) (see Fig. 1), which causes “freezing” of dynamics-a phenomenon called decoupling by 
“bang-bang” control Figure 2 illustrates this effect showing the evolution of Bloch vector (43) for different parameters 
of the model. 


B. Example 2: Two-level system with “transverse” coupling 

Here we examine a more involved case of directions transverse to the constant external magnetic field which appear 
in the NMR spin-bath coupling. The system’s Hamiltonian is the same as in the previous case (26), and the “magic 
angle” condition A = n/2 still applies. The Floquet operator U(T) and the evolution operator U(t) are also the same. 
The interaction Hamiltonian H int is now, however, of a form 

-Hint = (J 1 (D Ri + cr 2 (g) R 2 , (48) 

and in the context of NMR describes the “spin-lattice” coupling. 


1. The case of no control, f(t) = 0. 

In the absence of external control, the two terms in (48) lead to similar effects and act additively. For a thermal 
bath at the inverse temperature /3, TLS thermalizes with two decay times: Ti-for the diagonal elements (in a 3 basis) 
and 2Ti-for the off-diagonal elements of the density matrix. T\ is determined by the joint spectral density of a bath, 
satisfying the Kubo-Martin-Schwinger condition 

2 00 

7-l( w ) = X] J e~ lUlt (Rj(.t)Rj) B dt, 7j_(-w) = 


( 49 ) 


9 



FIG. 2. (Color online) Trajectories of TLS’s Bloch vector x(t) for longitudinal coupling and Lorentzian bath spectrum, (a) 
Example of evolution with parameters T 2 = 6.5 x 10 -3 , r c = 18.7, cjo = 18.8, and T — 1.5. (b) T 2 is reduced by a factor of five 
(T 2 = 1.3 x 10 3 ), producing a more unstable evolution and a faster approach to the stationary maximally mixed state, (c) 
The bath correlation time r c is multiplied by a factor of five (r c = 120) keeping T 2 constant, (d) The period of modulation T 
is decreased five times (T = 0.3)-both panels (c) and (d) result in much more stable evolution with long decoherence time. All 
plots were generated under A = 0 condition. 
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as 


Ti = [(l + e-*»)7j.(a;o)] 1 (50) 

The standard model for the spectral density of the acoustic phonon bath is given by the regularized expression 

( 51 ) 

where A is an effective coupling strength and uj cut is a cutoff corresponding to the Debye frequency. The similar 
formula can be used for the case of electromagnetic bath with dipole coupling which can be relevant for TLS models 
describing quantum dots, superconducting qubits, etc. 

If a bath with a spectral density of the type (51) is also coupled by the a 3 term; it does not contribute to the 
decoherence time, because 7_l(0) = 0. Therefore, under the joint influence of both types of baths the “spin-lattice” 
relaxation time is given by Eq. (50) and the decoherence time is finally given by 

^l = Y 2 + W 1 = \ ( 7 " (0) + (1 + e-/3 “°)^M). (52) 

2. The case of periodic kicks at resonance and zero temperature 

Now it is possible to discuss the case of periodic kicks in a similar fashion to the previous section, using the spectral 
density (51). However, formulas obtained for the relaxation times are very complicated and involve a series of special 
functions. Hence, we restrict ourselves to the resonance driving and zero temperature of the bath, i.e., A = 0, /3 = oo 
(Fig. 3). This case is sufficient to present the basic feature of the periodically kicked system-suppression of dissipation 
and decoherence for Tl cj cut . 


Y(w) [s' 1 ] 



FIG. 3. Example plot of spectral density function (51) for the zero temperature case. Parameters were chosen such that A = 1 
and cjcut — 5. 

The Heisenberg picture dynamics in the rotating frame of a 2 leads to the Fourier decomposition 

oo 

a 2 (t)= J2 (s+(g)e it ( f+1,!! )+S!(g)e it (-T +1!fl )) (53) 

q= — oo 


with matrices 


4 



0 0 \ 

IT35 V’ 



0 

0 


t)- 


(54) 
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given again in Floquet basis. Due to the condition A = 0, a 1 is a constant of motion and hence, because 7 _l( 0 ) = 0, 
it does not enter the generator which, in interaction picture, reads 


C P = Xb 1 - {~^ +qQ ) - ^{Sl(-q)S 2 _(q), p}^j 

+ XD- 1 - +<2^) (s%(q)pSl(-q) - (~q)Sl(q), p}^ , 


and in the explicit matrix form is the same as (36) 


Cp = -ru_(Q) 


P 11 ~ P 22 


-p 12 


P21 P22 ~ P 11 


with the single relaxation rate rj |_(D) given by the formula 


AO 3 


Q 


?7_l(D) = coth ( sinh 




47r 2 


2^cut \ 2 cj c 


(55) 


(56) 


(57) 


The plot of the relaxation rate (57) in Fig. 4 shows the suppression of decoherence for high kicking frequencies, 
D cj cu t. 

After returning to the original a 3 basis and to the Schrodinger picture (13), the full dynamics of p(t) is then 
expressed as [p = rj±(Q)] 


Pn(t) = \ [l + (-l)L^Je'^i(0) - P22(0))] , 


P 2 i(t) = -e iuJot 


- 2r,t (pi 2 ( 0 ) + P 2 i( 0 )) - (-l)L^Je-^(p 12 (0) - p 21 (0)) 


P22(t) = 1 ~ pll(t), Pl2(t) = p2l(t). 
Equivalently, the Bloch vector representation reads 


x\(t) = e 2r]t xi(0) cosuiot — ( — 1) L t J e ^^(0) sincj 0 t, 
X 2 (t) = e~ 2r]t xi(0) sincjo^ + (-1) LtJ e~ r]t X 2 { 0 ) cosc^o^ 
x 3 (t) = (-l)l±\e~* t x 3 (0). 


(58a) 

(58b) 

(58c) 


(59a) 

(59b) 

(59c) 


This evolution of Bloch vector coincides with Eq. (43) under the resonance condition A = 0, and again consists of 
three components: (a) rotation with the Larmor frequency uj o , (b) exponential decay with two decay rates, rj\\(T) 
and 2?7||(T), and (c) inversion of the slowly decaying component performed at the times nT. The exponential decay 
of all components of Bloch vector drives the TLS to the final maximally mixed state. 

Assuming that the parallel and perpendicular couplings are statistically independent their combined effect is de¬ 
scribed in the interaction picture by the simple generator of the form (36) with the overall decay rate rj = rj\\ + 


IV. CONCLUSIONS 

The formalism of Markovian master equations for periodically controlled open quantum systems weakly coupled to 
stationary environments allows us to revisit the well-known theory of “bang-bang” control providing a more detailed, 
quantitative and mathematically consistent description. Closed formulas can be applied directly to experiments 
concerning periodic control for various realizations of TLS. The methods are presented in a way which can be easily 
extended beyond the special choice of parameters used here to illustrate the main features of the pulsed control. 

In the case of NMR, by comparing the measured decoherence time for two values of kicking period: T r c and 
T < t c , one can determine the relaxation time of the bath, r c , by using Eq. (37). This is a new parameter beside the 
intensity of a signal and two relaxation times: “spin-lattice relaxation time”-Ti, and “spin-spin relaxation time”-X 2 . 
Such a new parameter, which depends on the environmental properties independent of those determining Xi and X 2 , 
could, in principle, increase the contrast in NMR imaging. 
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FIG. 4. Case of coupling by a 1 and a 2 with A » 0, A = 1 and uj CVlt = 1; comparison between functions rj±(Cl) (black) and 
7_l(Q) (gray, dashed). 
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Appendix A: Notes on Floquet formalism for kicked dynamics 


Here we sketch a more detailed, but not too formal derivation of the Floquet operator of a periodically kicked system 
as well as the evolution operator U(t), justifying formulas (16) and (18). In order to do so, we accept an approach 
mentioned in Ref. [25], namely, we consider a Floquet operator V (T) as an evolution operator, which propagates state 
V’(nT + ) of a system from a time just after the nth kick to a state fj((n + 1)T + ) of a system at time just after the 
kick n + 1. Here, the phrase “just after” is to be understood in a sense of limiting procedure, i.e., £ + will be defined 
as t + e for e being positive and arbitrarily small. Equivalently, one can choose another time frame for propagator 
with time 0 chosen “just after” the first kick, i.e., at t = e. In this new time frame, all kicks are slightly displaced, by 
e, in the direction of negative t. 

It must be noted, however, that many sources use seemingly different approach [22, 23, 25, 26] where U(T) is 
calculated as the evolution from state nT to ^((n + 1 )T _ ) with t~ denoting time just before time £, e.g., 
t~ = t — e with e again positive and small. This, being a rather opposite convention results in different Floquet 
operator, which turns out to be completely equivalent in most applications to ours. 

We consider evolution operator defined over interval [0 + ,£], such that the kick occuring at 0 is not counted. All 
subsequent kicks take place at times nT, n G Z. Informally speaking, each one of the 5 functions in the Dirac Comb 
XlfcL-oo ^ — fcT) can be approximated by a suitable well-behaved function, nonzero over the interval [nT - ,nT + ], 
with very high and very narrow peak at t = nT~ and being zero on the interval [nT + , (n + 1)T _ ] (such reasoning 
can be made rigorous by appropriate regularization of the S function). It is evident, that, in [0 + ,T _ ], the time- 
dependent part of Hamiltonian H(t) is not acting, so the system evolves freely via e ~ lHot . In the infinitesimally short 
interval [T _ ,T + ], where a kick takes place, the system experiences an infinitesimally short and infinitely strong kick 
by operator A A in such a way, that evolution is rapidly modified by an “instantaneous unitary” of the form 


exp 


T + 


-i j (H 0 + A AS(t - T)) dt 


= e lXA , 


(Al) 


where the equivalence comes from the fact, that during such short time only 6(t — T ) contributes significantly and Hq 
can be entirely omitted. Therefore, formally considering the limiting procedure one can replace t + and t~ with t and 
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the evolution over interval [0, T] is just a Floquet operator 

U(T) = e~ ixA e~ iH ° T . (A2) 

For graphical explanation, see Fig. 5. Now let us take t large, i.e., t can be expressed as t = nT + St for n G N and 



FIG. 5. Schematic interpretation of evolution on the interval [0 + ,t]. Each kick is taking place just before time nT and is 
included within evolution on [nT + , (n + 1 )T + ]. A whole evolution is then expressed as a composition U{t) = e ~ xH ^ 5t jJ {T) n 
of free evolution (generated by Ho) taking a time St with n copies of Floquet operator, which together give the remaining 
evolution over time nT. 

0 ^ St < T. Then, one has 


n = 


t 

T 


St = T 





(A3a) 

(A3b) 


with [x\ and {x} denoting the integer part (floor function) and fractional part (sawtooth function) of x with property 
{x} = x — [x\. The full evolution operator U(t) can therefore be composed of a free evolution over time St with n 
subsequent evolutions U(T) over time intervals [fcT, (k + 1)T], i.e. 


U(t) = e-iHoStjj^n 

(A4) 

which, after using Eq. (A3), is 


U(t) = e~ iHoT {^U{T)\-T\. 

(AS) 

Of course U(T) = e~ lHT and \_x\ = x — {x}, so one has 


U(t) = e~ l ** T lT \ 

(A6) 



_ e -iH 0 T{±} e iHT{±} e -iHt 



= P{t)e~ im , 

where operator P{t) = is clearly periodic: 

P{t + T)= e -«on¥} e ^{¥} (A7) 

= e -iH 0 T{j, + l} e iHT{j, + l} 

= e -iH 0 T{i} e iHT{i} ) 


since {x + 1} = {x} (see Fig. 6). 
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